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ABSTRACT 

A new mechanism for acceleration and enrichment of ■^He during impulsive solar flares is pre- 
sented. Low-frequency electromagnetic plasma waves excited by the Electron Firehose Instability 
' (EFI) can account for the acceleration of ions up to 1 MeV amu~^ energies as a single stage pro- 

cess. The EFI arises as a direct consequence of the free energy stored in a temperature anisotropy 
(T|j^ > Tji) of the bulk energized electron population during the acceleration process. In contrast 
to other mechanisms which require special plasma properties, the EFI is an intrinsic feature of 
the acceleration process of the bulk electrons. Being present as a side effect in the flaring plasma, 
these waves can account for the acceleration of ■^He and ''He while selectively enhancing ^He due 
to the spectral energy density built up from linear growth. Linearized kinetic theory, analytic 
models and test-particle simulations have been applied to investigate the ability of the waves to 
accelerate and fractionate. As waves grow in both directions parallel to the magnetic field, they 
can trap resonant ions and efficiently accelerate them to the highest energies. Plausible models 
have been found that can explain the observed energies, spectra and abundances of ^He and "^He. 

Subject headings: Sun: flares - Sun: particle emission - Sun: abundances 

1. INTRODUCTION McGuirc, & von Rosenvinge 1986). The divi- 

sion into these two categories can be done on the 
^ Solar flares are commonly divided into two ^asis of the duration of their soft X-ray emis- 

different classes: impulsive and gradual (Cane, 
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sion (Pallavicini, Scrio, & Vaiana 1977). But it 
is not only the timescale of the events that jus- 
tifies the distinction: the energetic particles ob- 
served in space from impulsive flares exhibit strong 
abundance enhancements over coronal values (Lin 
1987; Reames 1990, and references therein). Im- 
pulsive flares are usually dominated by energetic 
electrons and are characterized by ^He/''He ratios 
at 1 MeV amu~^ energies that are frequently 3 to 
4 orders of magnitude larger than the correspond- 
ing value in the solar corona and solar wind where 
^He/*He ~ 5 • 10"^. They also exhibit enhanced 
""^He/H and Fe/C ratios. Although the occurrence 
of "^He and ^^Fe enrichments are correlated in im- 
pulsive flares, the ratio ^He/Fe shows huge vari- 
ations as observed by Mason, Dwyer, & Mazur 
(2000). This suggests, that different mechanisms 
are responsible for the acceleration of the two 
species. Gradual flares usually have large ener- 
getic proton fluxes, small ^He/H and do not show 
large ^He/^He or Fe/C enhancements in the ener- 
getic particles, although approximately 5% admix- 
ture of suprathermal remnant particles from im- 
pulsive flares have been observed in gradual events 
by Tylka et al. (2001). The standard interpreta- 
tion for these observations is that the energetic 
particles in impulsive events origin in the energy 
release region on the sun while the energetic par- 
ticles in gradual events are accelerated via shocks, 
either coronal or interplanetary (Lin 1987; Luhn 
et al. 1987). 

Abundance ratios therefore are a valuable di- 
agnostics for the flaring plasma itself and in par- 
ticular the specific acceleration mechanism for the 
energetic particles. The selectivity of the mech- 
anism, especially for ^He and ^He indicates res- 
onant processes such as gyroresonant interaction 
of plasma waves with the ions. Theoretical ideas 
therefore focus on the unique charge-to-mass ratio 
of ^He which allows it to be selectively pre-heated 
or accelerated via gyroresonance. 

A well-known theory of the initial set among 
theories for "'He enhancement was published by 
Fisk (1978), explaining the preferential accelera- 
tion of the ions by electrostatic ion cyclotron (EIC) 
waves at a frequency in the vicinity of the gyro- 
frequency of '^He . The waves are excited by an 
electron current and interact with ^He via cyclo- 
tron resonance. A large enhancement of ^He/H is 
required in the ambient plasma for this instability 



to excite waves above the ^He gyrofrequency. 

More recently a theory was suggested by 
Temerin & Roth (1992) who accounted for the 
preferential "^He acceleration proton electromag- 
netic ion cyclotron (H+ EMIC) waves. These 
waves are driven unstable by non-relativistic (keV 
range) electron beams and their frequencies lie at 
around the "^He gyrofrequency at almost perpen- 
dicular propagation. In Temerin & Roth (1992), 
auroral observations of keV electron beams and 
H"*" EMIC waves are taken as experimental evi- 
dence that 11+ EMIC waves also may acquire a 
substantial fraction (order of few percent) of the 
electron beam energy under coronal conditions. 
At the Sun, plasma emission from tens of keV 
electron beams on open magnetic field lines is the 
explanation for type III radio emission with its 
equivalent, the U-bursts, on closed field lines (in 
loops) at beam energies of order keV. In order to 
excite H+ EMIC waves fulfilling the requirements 
of this model, electron beams of much higher den- 
sity than the observed ones have to be postulated 
(Miller & Vinas 1993) with no direct observational 
evidence. Moreover, it is difficult to explain from 
a theoretical point of view that the free energy in 
the electron beam is transferred to the H+ EMIC 
waves and not to the much faster growing (^ 4 or- 
ders of magnitude in growth rate) electron plasma 
waves. 

In this work an alternative model for accel- 
eration of ■'He and its enhancement over ^He is 
presented. The approach is different from the 
models described above. While the models men- 
tioned above postulate rather special plasma prop- 
erties (^He/H overabundance in the pre-flaring 
plasma, dense low-energy beams) in order to pro- 
duce the required plasma waves, the model pre- 
sented herein explains the unique overabundance 
of ^He by plasma waves excited as an intrinsic 
feature of the electron acceleration process itself. 
Parallel propagating, lefthand polarized electro- 
magnetic waves driven by the Electron Firehose 
instability (EFI) can account for ''He accelera- 
tion via gyroresonant interaction. As suggested 
in Paesold & Benz (1999) such Electron Firehose 
(EF) waves are excitcid by anisotropic electron dis- 
tribution functions (Ty > T±) that occur in the 
course of the acceleration process of bulk electrons 
in solar flares. 

Although the EF waves are not narrowbanded 
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around the gyrofrcqucncy of "^Hc as the H+ EMIC 
proposed by Temerin & Roth (1992), selectivity 
of the process is achieved by the natural profile of 
the spectral wave energy. No additional assump- 
tions besides the electron anisotropy are needed, 
and the model can be embedded as an intrinsic fea- 
ture in acceleration scenarios such as transit-time 
damping, the currently most popular stochastic 
acceleration model. 

In § 2 the basics of the new acceleration scenario 
are described. The properties of the EF waves un- 
der coronal conditions are presented in § 3.1. An- 
alytical and numerical results on the heating rates 
are presented in § 3.2 and § 3.3. The mechanism 
for enhancement of ^He over ^He is described in 
§ 4, and heavier ions are discussed in § 5. § 6 
concludes this work. 

2. BASIC IDEA 

Among the most promising scenarios for accel- 
erating electrons to observed energies in impul- 
sive solar flares is transit-time damping accelera- 
tion (Fisk 1976; Stix 1992), the magnetic analogon 
of Landau damping. The following scenario was 
presented by Miller et al. (1996): Electrons are 
accelerated from thermal to relativistic energies 
by resonance with low-amplitude fast-mode waves 
having a continuous broadband spectrum. The 
magnetic moment of the particle interacts with 
the parallel gradient of the magnetic field as in the 
well known Fermi acceleration (Fermi 1949; Davis 
1956). Contrary to the classic Fermi process, 
transit-time damping involves small-amplitude 
magnetic compressions. While Fermi acceleration 
is the result of large numbers of particles being 
reflected by randomly moving magnetic compres- 
sions, transit-time damping is a process of rather 
resonant nature. In the limit of very small am- 
plitudes, only particles of the same parallel speed 
as the wave phase velocity can be reflected by the 
magnetic compressions, i.e. = to/k which is 
the Landau resonance condition. Interaction with 
a wave changes the particle's parallel speed and 
therefore allows it to interact with another wave 
of the continuous spectriun. In the average this 
process results in stochastic acceleration. 

The wave spectrum is assumed to origin from 
cascading fast-mode waves, initially excited at 
very long wavelengths as a direct output of the 



primary energy release in impulsive flares. The 
cascade channels the released energy through an 
inertial region to fc- values small enough to acceler- 
ate electrons out of the thermal background pop- 
ulation. 

Other acceleration scenarios operating in im- 
pulsive solar flares have been proposed and can 
roughly be divided into three categories: Shock 
acceleration, acceleration by parallel electric fields 
and stochastic acceleration by MHD turbulence 
including the model described above. A detailed 
review can be found in Miller et al. (1997) and 
references therein. 

Even though the nature of the actual accelera- 
tion mechanisms is unclear, all the above mecha- 
nisms that can account for accelerating the bulk of 
electrons up to the observed energies of ^ 20 keV 
have an important feature in common: Particles 
are preferentially accelerated in parallel direction 
with respect to the background magnetic field. 
Parallel dc electric fields trivially accelerate only 
along the magnetic field while stochastic scenar- 
ios as transit-time damping act via small ampli- 
tude magnetic mirroring which is only capable of 
transferring energy in parallel direction if no addi- 
tional scattering mechanism is provided (Lenters 
& Miller 1998). Works of e.g. Wu (1984) and 
Leroy & Mangeney (1984) describe the parallely 
directed energization of electrons at the earths 
bow shock via shock drift acceleration with quasi 
perpendicular shocks. 

The distribution function of the accelerated 
particle species in velocity space therefore is ex- 
pected to become more and more anisotropic in 
the course of the energization process. While the 
perpendicular temperature remains virtually con- 
stant during the acceleration, the parallel velocity 
of the particles increases. If the energization in 
parallel direction is from a thermal level of some 
0.1 keV up to 20 keV or more, the anisotropy 
in parallel direction is substantial. The plasma 
therefore is modeled by a bi-maxwellian with dif- 
ferent temperatures in parallel and perpendicular 
direction. In doing so it is assumed that no effi- 
cient scattering mechanism is present that could 
maintain isotropy on acceleration timescales. Al- 
though it is not expected that acceleration retains 
a bi-maxwellian distribution function, it is a good 
approximation in order to describe propagation 
properties of plasma waves. As shown in Paesold 
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& Bcnz (1999) the EFI can occur in such a situ- 
ation and give rise to EF waves. The parallel EF 
waves are purely transverse electromagnetic and 
lefthand circularly polarized. The sense of polar- 
ization and can change in certain k ranges if a 
similar anisotropy of the protons is assumed. In 
the following the EF waves arc lefthand polarized 
if not mentioned otherwise. A representative dis- 
persion plot is depicted in Figure 1. 

Taking into account the uncertainties in the ac- 
celeration region, including possible pre-heating 
mechanisms, reasonable pre-flaring plasma con- 
ditions of an impulsive flare range within Bq « 
100 — 500 G for the background magnetic field, 
Ua « 10® — 10^^ cm~^ in number density and 
about Ta ~ 10*' — 10^ K for the proton and elec- 
tron temperature (Pallavicini et al. 1977). For 
the numerical example herein the following pre- 
flaring plasma parameters are chosen: = n.p = 
5-10^° cm-^T| II = II = 1-10^ K,Bo = 100 G. 

Under these exemplary conditions the EF waves 
propagate below about 3 • Oh- Several ion gyro- 
frequencies are indicated in Figure 1 and they lie 
well within the wave spectrum. Resonant acceler- 
ation occurs if the condition w — k\\V\\ — Zf2/7 = 
is satisfied. Here, W|| and 7 are the parallel parti- 
cle speed and Lorentz factor, and Q. is the gyro- 
frequency of the according particle. This condition 
is well satisfied for ions like 3Hc2+ and '^}lc^+ . The 
model presented herein assumes acceleration of the 
ions via the most effective gyroresonance at Z = 1, 
the so called cyclotron resonance. 

Due to the symmetry of the distribution func- 
tion with respect to uy, the dispersion of the 
EF waves is the same for the transition from 
k — » — k while keeping lefthand polarization and 
a; > (Hollweg & Volk 1970). The resulting 
wavefield from the EFI therefore consists of waves 
propagating parallel to the magnetic field (pos- 
itive k-branch) as well as anti-parallel (negative 
k-branch). This property of the EF waves is of 
special interest and crucial for efficient acceler- 
ation: A unidirectional wavefield of electromag- 
netic waves always exhibits a force parallel to the 
background magnetic field pushing the particle out 
of resonance, limiting the acceleration time and, 
hence, the reachable energies. While other mod- 
els (e.g. Roth & Temerin 1997) solve the problem 
by imposing a background magnetic field geome- 
try, i.e. a field gradient, to force the particle to 



regain resonance, this is not necessary for the case 
of the EF waves: resonant acceleration by counter- 
propagating electromagnetic waves leads to oscil- 
lating parallel forces and does not drive the par- 
ticle permanently out of resonance. In contrary 
they naturally 'trap' the particle in the wavefields 
by bouncing it from resonance with the positive 
branch to resonance with the negative branch of 
the wavefield and vice versa. This ping-pong trap- 
ping strongly enhances the total resonant interac- 
tion time of the particle with the wave and allows 
significant acceleration in perpendicular velocity. 
The mechanism works even for small wave fields 
where 'normal' wave-particle trapping is not effi- 
cient. Similar situations have been investigated by 
e.g. inferring counterpropagating Alfvcn waves to 
the problem of solar flare proton acceleration (Bar- 
bosa 1979) and cosmic-ray acceleration (Skilling 
1975). This acceleration can be regarded as a spe- 
cial case of wave turbulent stochastic acceleration 
with waves explicitly counterpropagating in one 
dimension. 

The differences in enrichment of '^He and ^He is 
the result of varying growth of the waves at differ- 
ent frequencies. As the waves grow from a thermal 
level to a saturated state, the spectral energy den- 
sity develops according to the growth rates from 
linear theory depicted in Figure 1 . Since the waves 
are excited nonresonantly by the electrons, the 
saturated wave energy spectrum can be approx- 
imated from the profile of the linear growth rate. 
This is not true for resonantly excited instabilities, 
where wave modes can still grow while others have 
already saturated. If each mode is in resonance 
with only the particles fulfilling the resonance con- 
dition, some modes can grow longer than others 
and the saturated energy spectrum of the waves 
has to be determined by other methods. How- 
ever, if the instability is nonresonant as the EFI, 
each mode grows according to the total free en- 
ergy available (the driver is basically a pressure 
anisotropy of the electrons) and the whole elec- 
tron population contributes to the growth of each 
mode. All modes therefore saturate at the same 
time, namely when the pressure anisotropy drops 
below the threshold, and the spectral energy is 
frozen at that time. The wave dispersion is al- 
tered by the decreasing anisotropy as result of the 
erosion of the particle distribution. However, a 
former analysis shows (Paesold & Benz 1999) that 
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this docs not severely alter the slope in the growth 
rate around '^He an ''He. 

Although the differences in growth rates are 
quite small in the region of interest (see Fig. 2), 
the differences in wave energy will become signif- 
icant after several growth times r = I/7. As a 
rough estimate one obtains, by assuming constant 
growth rates in time, that two modes with ini- 
tial growth rates 7a = 7Max and 7b = 0.9 7Max 
and equal energy content have a ratio in energy 
of Wb/Wa ~ 0.1 after a time r = 7^^ of only 
Ri 10. This estimate illustrates the rather strong 
influence of growth rate on the resulting wave en- 
ergy level and, hence, on the selectivity of cyclo- 
tron acceleration due to the 7 profile in wave fre- 
quency. The analysis in Paesold & Benz (1999) 
shows that the positive slope in the growth rate 
profile around the ^He and *He gyrofrequency that 
is needed to selectively enhance '^Hc is a stable fea- 
ture of the EFI and not very sensitive to changes 
in the plasma parameters for solar pre-flaring con- 
ditions. 

In the following it will be established that the 
mechanism meets the following requirements for 
the observed ^He enrichment: 

(i) The accelerator has to energize the ions 
above ~ 1 MeV amu~^ on timescales of 1 s. 

(a) ■^He/'^He > 0.1 above energies of about 
1 MeV amu~^ . 

(in) A time integrated total of about 10'^'^ "^He 
nuclei is required to account for the observed par- 
ticle fluxes (Reames et al. 1994). 

(iv) The ^Hc spectrum is harder than the ^He 
spectrum in the range of 0.4 — 4.0 MeV amu~^ 
(Moebius et al. 1982). 

(v) ^He exhibits a turnover in the particle en- 
ergy spectrum at around a few ~ 100 keV amu~^ 
(Mason et al. 2000). 

3. ACCELERATION 

3.1. Electron Firehose Wave Properties 

Parker (1958) pointed out that a magnetized 
plasma with a pressure anisotropy in parallel di- 
rection to the magnetic field can become unsta- 
ble to low frequency Alfvcn waves. This instabil- 
ity is known as the Firehose instability and is of 
completely non-resonant nature: neither the elec- 
trons nor the protons are in resonance with the 
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Fig. 1. — Dispersion relation of Electron Fire- 
hose waves calculated from linear theory. The fre- 
quency Ur (solid curve) and growth rate 7 (dashed 
curve) of the EF waves are shown normalized 
to the proton gyrofreqiiency as function of par- 
allel wavenumber A;|| normalized to the inertial 
length of the protons. The plasma parameters are 
T|/rX = 15, ^/Tl = 2, Tl = Tl = 1 ■ 10^ K, 
Ue = Up = 5 - cm-3 and Bq = 100 G. The 
sign convention is that fc > and a; > mean left- 
hand circularly polarized; 7 > refers to growing 
modes. The dotted horizontal lines indicate the 
gyrofrequency of the according ion species. 

excited waves. An extension of this instability 
to higher frequencies was presented by HoUweg 
& Volk (1970) and Pilipp & Volk (1971). This 
branch of the instability has been termed the Elec- 
tron Firehose Instability. Here, the bulk of the 
protons is resonant with the waves which are non- 
rcsonantly excited by the electrons. The electrons 
are anisotropic and drive the waves while the pro- 
tons carry the wave. A typical dispersion relation 
of the EF waves is displayed in Figure 1. The 
dispersion of the EF waves is computed by IDL- 
Whamp (Paesold 2002), an easy to use IDL in- 
terface to the WHAMP code originally developed 
by Ronnmark (1982). The code provides the user 
with the full solution of the dispersion equation in 
linearized kinetic theory. 

The influence of the presence of other major- 
ity ions such as *He on the EF wave dispersion 
is negligible. For a plasma consisting of 5% ^He 
and 95% H the wavenumber of maximum growth 
and the frequency at maximum growth is shifted 
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by values of order of 1% with respect to the val- 
ues of a pure H plasma. Ions other than H+ have 
therefore been omitted in the following when com- 
puting the dispersion relation of the EF waves. 
Due to the assumed anisotropy of the protons 
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Fig. 2. — Growth rate vs. frequency w^/fl/f 

for the same plasma as in Figure 1. Vertical dotted 
lines indicate the gyrofrequencies of the according 
ions. 
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Fig. 3. — Maximum growth rate 7Max/f^ff vs. fre- 
quency ijj/Q.H- For each of the three lines the 
electron anisotropy is constant and the proton 
anisotropy is varied. The most upper line (tri- 
angles) depicts r|[/T| = 20, the middle Unc (di- 
amonds) T^/T1 = 15 and the bottom line (aster- 
isks) T^\/Tl_ = 14. The small numbers refer to the 
according values of /T^. The vertical dotted 
line indicates the ^He^"*" gyrofrequency. 



the mode is righthand polarized {uir < 0) at small 
values of k in Figure 1. According to HoUweg & 
Volk (1970) the condition for the change of po- 
larization sense to occur is an additional proton 
anisotropy /T^ ^ 2. Such a proton anisotropy 
could result e.g. from transit-time acceleration 
as described for the electrons in § 2. Note, that 
the proton anisotropy is not a necessary condi- 
tion for instability. It has been taken into account 
herein only for the sake of generality. The de- 
pendence of the maximum growth rates on pro- 
ton and electron anisotropy can be seen in Fig- 
ure 3. The proton anisotropy has been varied from 
1 to larger values for three values of the electron 
anisotropy. Only lefthand modes have been taken 
into account. With increasing proton anisotropy 
the frequency at maximum growth approaches and 
finally crosses the ^He gyrofrequency. With de- 
creasing electron anisotropy less anisotropic pro- 
tons arc needed to shift the most growing mode to 
the ^He gyrofrequency. This shows that for rea- 
sonable anisotropies the slope in the growth rate 
profile in the vicinity of the gyrofrequency of "'He 
is always positive, such that waves around the '^He 
resonance grow faster than waves around the *He 
gyrofrequency. This is a very stable feature of 



the EF wave spectrum and can be expected, once 
above the instability threshold, for a large range of 
plasma parameters that allow instability and rea- 
sonable anisotropies of protons and electrons (for 
a detailed analysis of EFI thresholds and spectra 
in dependence on plasma parameters see Paesold 
& Benz 1999). 

An additional mode can be found at oblique 
propagation angles (Paesold & Benz 1999; Li & 
Habbal 2000). For the reported examples, the 
oblique mode grows faster than the parallel mode. 
In addition to purely growing waves, the oblique 
mode contains waves that may resonate with coro- 
nal ions. Relevant for the flare application and 
mode comparison is the instability threshold and 
the non-linear evolution. Also unknown in this 
context are the effects of ion anisotropy, power- 
law electron distribution and other parameters. 
Thus we concentrate here on the well known par- 
allel mode and leave the exploration of the oblique 
mode to future work. 
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3.2. Heating Rates 

3.2.1. Unidirectional Propagation 

The ion's energy change is dominated by the 
resonant interaction with the wave's electric field, 
which is decomposed into its parallel component 
£^11 = E^iicos^*, where \E' = k\\z + k±x — tot 
is the wave's phase at the location of the ion, 
and two circularly polarized components E± = 
i?±(cos5',=Fsin5') with amplitudes E± = {E^ ± 
Ey)/2. In the Fourier domain we decompose the 
wave into a set of monochromatic waves sepa- 
rated by a frcqucnc;y interval Acj. In one of these 
monochromatic waves with amplitudes ey and e± 
and frequency oj, the ion gains or loses energy ac- 
cording to 

W = qv_\_e±cos4>±-\- qv\\e\\cos'^ , (1) 

where the relative angle (j)± = 6* ± \1/ between the 
wave electric field and the perpendicular velocity 
vector {vx,Vy) = w_L(cos6',sin0) of the ion, moving 
with the instantaneous gyrophase 9{t). 

The ion may be driven into or out of cyclo- 
tron resonance due to changes in the ion's paral- 
lel speed, which in turn depends on the wave's 
electric and magnetic field. The latter is ob- 
tained via Faraday's law B = —V x E as, (l/w) 
[±(fc||e±) sin^*, (fc||e±—A;_Le||) COS'S', T{ki.e±) sin^f], 
so that the parallel acceleration is 



mw|| = q 
-I- 



e± cos (f)±+ 



k±v± cos6 , „ 
1 1 eii cosM/ 



• (2) 



A possible mirror force due to a background mag- 
netic field gradient has been omitted since Bq is 
assumed to be homogeneous herein. 

The instantaneous ion angular frequency, 9, is 
given by the ratio of the force on the ion, perpen- 
dicular to the velocity and the background mag- 
netic field, and the moment perpendicular to the 
magnetic field. Considering only one monochro- 
matic wave at frequency u), the phase equation 
^± = 6±^ is given by 



-fi j ± (fc|| -y II — LO + k±v± cos < 

qe± ( k\\v\\ 

1 — I sm , (3) 



mv\_ \ LO 
where fi, is the gyrofrequency of the ion i. 



In case of the EF waves, where k± = 0, ey = 0, 
and e+ = 0, Eqs. 1 and 2 simplify to 



W 



qe-V± cos( 



qe- 

v± = cos ( 

m 



qe- kHV± 
v» = — cos ( 



(4) 
(5) 

(6) 



For a particle i in resonance with a monochro- 
matic EF wave the frequency mismatch parameter 
^ = rij -I- fc||i;|| — u! vanishes. The righthand side of 
Equation 3 therefore reduces to a nonlinear pendu- 
lum equation for the phase 0_ . In resonance </>_ {t) 
can have a stable solution when it stays close to 
0. In this case the particle is accelerated in par- 
allel direction until it is expelled from resonance 
(Eq. 6). For particles not in resonance, (j)- changes 
rapidly and no energy is gained in the time aver- 
age. 

In a set of N monochromatic waves, numbered 
by i = 1,...-/V with equal amplitudes e_j- = 
-Erms/ and equally spaced by a frequency Aiv 
with Scj = NAu), the ion's perpendicular motion 

has a phase (j> = — with respect to each mode. 
This gives a set of N phase equations 



Pj = -rii - kiijvi\ - Uj 



^ TP 

sr^ qE„ 



smc 



(7) 



and Equations 5, 6 have to be summed over the 
contributions of the N modes. 

The equations of motion are in general not in- 
tegrable because of the sin (j)i terms in Equation 7. 
However, a maximum heating rate and the typical 
time after which an ion is driven out of resonance 
by the mean parallel force of the uni-directional 
propagating waves are estimated. The latter lim- 
its the total energy gain of the ion. 

Within a time interval t, the energy gain of the 
ion is mainly due to the wave modes in a band- 
width 2'k/t around the frequency « 17^ -f fc||W||. 
The wave may be described by a mode spacing 
Aw ~ 2-k/t, yielding e_ k, Ermsl \/ Sojt /2t:. At 
best, the ion stays in resonance with one mode jo 



such that 



over the whole time period r. 



This gives, using Equation 5, 

V^{t + T)-V^{t) 



(8) 
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and if t is sufficiently large so that v±{t) <C v±{t+ 
t), the heating rate per mass is at most 



VJ_VJ_ 

m 



2 

rms 



Mev 



6.1 



amu s 



(9) 

(10) 



100 V 



lO^rad s" 

5uj 



A 



independent of the length of the time interval r. 
Thus, Equation 9 predicts a linear increase of en- 
ergy in time. The ratio Q/A is the charge-to-mass 
ratio of the ion in atomic units. 

The resonant heating is limited to the time, 
during which the ion is in resonance with one of 
the wave modes. The ion dynamics are dominated 
by the timescales for parallel and perpendicular 
acceleration 

v\\ rnvnio 



mv± 



(11) 
(12) 



The initial conditions are usually ^;_lo ~ v\\o 
with fc||f|| <C Su) and r_Lo/T||o < 1. With increasing 
v±, the timescale t± increases and the timescale 
Ty decreases. At time t ^ tr, the ion eventu- 
ally reaches a parallel speed |t;||| of order 5w/k\\ 
and is driven out of resonance. The acceleration 
timescale at tr is 



U)T» 



2TTqErmskf,V±{tr) 



(13) 



Estimating tr shows that the total time of the res- 
onant acceleration process is not much longer than 
Ty^r, and thus from Equation 13 follows 



.3 /. -X 2 



1 6— (—\ (-\ (^ lOOV "1" '' ' 

LV \ CO J \QJ \ Erms 
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With Equation 9 the corresponding energy gain 
is v'^ ~ w^w/Zcy, which is only a few keV amu~"^ 
for the given parameters. It is too little to explain 
the observation, assuming realistic conditions. 

3.2.2. Bidirectional Propagation 

Further acceleration is possible in the situation 

of the EFT where counterpropagating electromag- 
netic waves of the same sense of polarization are 
present. The basic idea is, that the time average 
of the wave forces in parallel direction of counter- 
propagating waves averages out to some extend. 
That way, the ion can stay in resonance for much 
longer time. 

In the following lefthand circular polarization is 
assumed. Plus and minus signs in index positions 
therefore indicate the direction of propagation and 
no longer sense of polarization. The parallel accel- 
eration now reads 

N 



^ej cos$+^j— Cj cos$_j- 



(15) 



where the phases = 9 T ^llj^^n* + <^jt are 

introduced. Than can be expressed by 
and we rewrite Equation 15 as 



N 



qk\\jv± 



{ej - ej)cos$+„ 



+ 2 Bj sin [6 + Wjt) sin (fc|| . (16) 

The first term in the bracket is similar to the uni- 
directional case (Eq. 6) except that the difference 
between the bidirectional fields enters now and re- 
duces the parallel acceleration. As in the unidi- 
rectional case this term is nearly constant in res- 
onance, i.e. Ri 0. The second term in the 
bracket contains a product of two sinus terms. At 
times ty < 27r/fc||^jU|| there is a mean change in v\\ 
when averaging over t due to the first sinus term. 
The second sinus term, sin (A;|| leads to an 
oscillatory motion in for ty > 27r/fc||jW|| and 
thus slows down the escape from the resonance 
region and, eventually, allows the regain of reso- 
nance. The maximum energy, that can be reached 
by the ion therefore can be increased. The dynam- 
ics of U|| (Eq. 16) is nonlinear and treated in nu- 
merical simulations with the exact wave spectrum 
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and dispersion relation of the counterpropagating 
waves. 

3.3. Test-Particle Simulation 

The rclativistic equations of motion for a par- 
ticle of charge q and mass to in a field of N waves 
and a homogeneous background magnetic field Bq 
are given by 

, N 

^=q^xBo + qJ2(E, + ^xB,) (17) 
fe=i 



where x is the particle position vector, v the ve- 
locity and Ek and B^ are the electric and mag- 
netic fields of the wave k. 

Test-Particle trajectories are calculated by in- 
tegrating a dimensionless form of Equations (17) 
and (18) with a standard leap-frog mover following 
Birdsall & Langdon (1991) (half acceleration with 
an extrapolated electric field, rotation around the 
instantaneous magnetic field, which includes the 
wavefield, and half acceleration with updated ve- 
locities). The code used herein is a leap-frog ver- 
sion of the code used by Miller & Viiias (1993) 
and has been tested against their version. The 
simulation results of both codes are in excellent 
agreement on timescales relevant for the analysis 
herein. 

The time step in the simulation was r = 
O.lfijj^. The background magnetic field is con- 
sidered to be homogeneous. The wave frequen- 
cies, parallel wavenumbers and relative strength 
of electric field components are obtained from the 
IDLWhamp code described in § 3.1. The electric 
field strength Ek of one single wave k is obtained 
from growth according to linear theory and the 
magnetic field components are calculated accord- 
ing to Faraday's law. A wavefield consisting of 
2000 monochromatic waves has been applied, con- 
fined to a range of frequencies where the growth 
rate does not drop below 60 % of the maxi- 
mum growth rate. All properties of the wavefield 
besides the relative phases are symmetric with 
respect to the zero of fcy. The phases were cho- 
sen randomly for all waves independent of the 
direction of propagation. Only waves propagat- 
ing parallel to the background magnetic field are 
considered. Some properties of the applied wave 



spectrum arc displayed in Figure 4. The waves 
are purely transverse and lefthand circularly po- 
larized. 

The spectral energy density of the wavefield is 
obtained by letting the waves exponentially grow 
according to W{t, k) = Wo{k) exp {2"t{k) ■ t) from 
a thermal level M^o(^) = up to an energy den- 
sity of W{t, k) at a time r and ^{k) is assumed to 
be constant during this time. The time of growth 
has been chosen such, that the total energy den- 
sity in the wave field is Wtot(T)/f/o = 0.01 at a 
time r, where Uq = Bq/Stt is the energy density 
of the background magnetic field. 

The plasma parameters are Ue = Up = 5 ■ 
10^° cm-3,rX = Tl = 1- 10"^ K, TpTl = 15, 
T||' /T^ = 2 (see Fig. 1) and the background mag- 
netic field is Bo = 100 G. 

To reduce computing time, the simulations have 
been divided in two parts: I) For each of the two 
ion species, a small number of particles have been 
followed for a long time (^ 1 s) to establish the 
ability of the process to accelerate ions to the re- 
quested energies (1 MeV amu~^ ). About 500 par- 
ticles were computed in order to determine also 
the spread in energy of the particle population. 
11) A larger number of ions (5000 particles for 
each species) have been followed for shorter times 
(~ 0.052 s) to establish the energy distribution 
which then was applied to the results of the first 
runs. 

A representative run of the part I) is depicted 
in Figure 5. The initial energy of the ion is 
100 eV amu""'^, corresponding to a thermal veloc- 
ity at a temperature of 1 • lO*" K, and the pitch 
angle cosine is /x = 0.5. The heating rates in the 
simulation are not very sensitive to variations in 
the initial energy of the ion. All particles start at 
space coordinates equal to zero. For each parti- 
cle the phases of the wavefield have been chosen 
randomly. The total energy in the wavefield is 
Uw/Uo ~ 0.01 corresponding to the wave electric 
field amplitudes Ek/Bo depicted in Figure 4. The 
plasma parameters are chosen as described in § 2. 

The typical energy evolution of a ^He ion shows 

intervals of zero gain (Fig. 5a). These arc times 
when the ion looses resonance with the wavefield 
due to a large excursion of from its oscilla- 
tory behavior described by Equation 16. Dur- 
ing these intervals, the ion oscillates in the wave- 
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Fig. 4. — Properties of the wave spectrum assumed for the test particle simulations, (a) Dispersion of the 
wavefield used for the test-particle simulation obtained from linearized theory. Each branch (positive and 
negative k) consists of 1000 monochromatic waves. The dotted lines indicate the resonance frequencies of 
the according ions, (h) The electric field amplitudes of the wave spectrum at the end of the growth time. 
The vertical dotted (dashed) lines refer to the k values corresponding to the according ion gyrofrequencies. 
(c) The corresponding spectral energy density. 



field without being significantly accelerated. Due 

to non-resonant interaction with the waves, the 
ion is scattered back into resonance and can be 
further accelerated. The loss and gain of reso- 
nance is mirrored in panel (h) and (d) of Fig- 
ure 5. During the intervals of no resonance (e.g. 
(t X I^h) X 10""^ « 8 — 15), the ion propagates freely 
and the coordinate along (Fig. 5fej is just a 
straight line, increasing or decreasing proportional 
with time. As can be seen in Figure 5d the paral- 
lel momentum no longer oscillates around zero in 
the above time interval. When the ion re-enters 
resonance, the parallel momentum starts oscillat- 
ing around zero, indicating ping-pong trapping of 
the particle, and the coordinate along Bo does not 
linearly change anymore. 

To further illustrate the ping-pong behavior, 
the so-called frequency mismatch parameter ^ = 



ujr — k\\V\\ — fis (where fis is the ^He^'*' gyro- 
frequency) has been plotted in Figure 6. Reso- 
nant interaction of the particle with a wave means 
^ 0, reproducing the resonance condition. Panel 
(a) of Figure 6 displays a small portion of the time- 
serie of Figure 5 a and the three vertical lines in- 
dicate times at which ^ was computed. As can be 
seen in Figure 66, the particle is out of resonance 
at times ti and ^3 while it is resonantly interact- 
ing with both wave branches at time t2- Accord- 
ingly, the parallel normalized momentum p^/'mc 
does not exhibit any oscillatory behavior around 
zero at times ti and in Figure 6a but clearly 
oscillates around zero at 

Acceleration is clearly very efficient with an 
averaged systematic energy gain (dE/dt) of « 
1.2 McV amu"-'^ s^^ (Fig. 5a). A linear energy 
increase is predicted by Equation 9. The paths of 
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Fig. 5. — Temporal history of a typical ^He^"*" test-particle, (a) The kinetic energy in kcV amu~^ of an 
initially thermal ^He ion (100 eV amu^^) vs. time in units of lO'^fl^^. (b) Distance traveled by the particle 
along the background magnetic field Bq. (c) Perpendicular and (d) parallel momentum of the particle vs. 
time. 



more than 500 ''He ions have been simulated and 
exhibit similar energization. 

4. ABUNDANCE ENHANCEMENTS 

Among the overabundant ions observed during 
impulsive solar flares, '^He takes a special posi- 
tion. Whereas most ions exhibit an enhancement 
by factors of 1-10 with respect to coronal values, 
•^He shows an excess of ~ 2000. Since it has been 
established that the selective enhancement of ion 
abundances correlates with the charge-to-mass ra- 
tio (Reames 1998) and hence the gyrofrequency 
of the according ions, resonant wave-particle in- 
teraction is generally believed to be the acceler- 
ating mechanism. This led earlier theories to the 
approach of finding special plasma waves existing 
only in a narrow range in the vicinity of the '^He 
gyrofrequency. In order to generate the appropri- 
ate waves (e.g. EIC waves (Fisk 1978), H+ EMIC 
waves in Temerin & Roth (1992)) peculiar plasma 



conditions have to be postulated. In the first case 
a pre-flare enhancement of ^Hc over H and in the 
latter case dense low-energy beams have to be as- 
sumed in order to excite the waves in the appro- 
priate frequency ranges. 

A different approach is investigated in the 
work presented herein. EF waves, excited by an 
anisotropic velocity distribution function result- 
ing from bulk acceleration of electrons, accelerate 
■^He and ■*He . Selectivity is achieved by an en- 
ergy density profile of the waves self-consistently 
produced from linear growth. For nonresonant 
instabilities, this profile is approximately true not 
only during the growth phase but also for the sat- 
urated state of the instability (see § 2). As can 
be seen from Figure 4 the rcsiilting energy density 
profile clearly supports an enhanced acceleration 
of ^He above ^He. The energy density in the wave 
field is higher at the gyrofrequency of "^He than 
*He. More energy is therefore available for the ac- 
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Fig. 6. — Zoomed region of tiie timeserie displayed 
in Figure 5 a. Panel (a) shows the time interval 
8.58-8.67 in normalized time of the parallel mo- 
mentum, (b) Plots of the frequency mismatch pa- 
rameter ^ = u)r — — O3 at three different times 
as indicated in panel (a). 

celeration of ^He. The quantitative heating rates 
of an ensemble of ions of the two species have 
been determined by numerical simulations where 
the total energy density in the wavefield has been 
treated as a free parameter. 

In order to compare the simulation results to 
the observational values of '^He/^Hc , abundances 
have to be compared at the same energies. Obser- 
vations from the ISEE 3 spacecraft were taken in 
the 1.3 — 1.6 MeV amu~^ channel (Reames et al. 
1994) and observations by Moebius et al. (1982) 
with the ISEE 1 and ISEE 3 spacecrafts were 
taken in the 0.4 — 4.0 MeV amu~^ range. Accord- 
ing to these measurements, a viable mechanism for 



•^He acceleration and its enhancement above ''He 
has to meet the requirements of § 2: (i) Ion en- 
ergies ~ 1 MeV amu"'' on timescales of 1 s, (ii) 
^Hc/'^Hc > 0.1 above 1 MeV amu"^ , (in) a to- 
tal of about 10"^^ ^He nuclei above 1 MeV amu~^, 
(iv) ^He spectrum harder than the ^He spectrum 
in the range of 0.4 — 4.0 MeV amu^"' and (v) a 
turnover in the particle energy spectrum of ^He at 
around 100 keV amu~^. 

In the following, the presented model is inves- 
tigated in view of these requirements. 




40 60 
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Fig. 7. — Temporal evolution of thermal energy 

Eth{t) as defined in Equation 19. A maxwellian 
distribution has been fitted to an ensemble of 500 
particles for each species. 

The ability of the EF waves to accelerate the 
ions to the relevant energies has been shown in § 3 
and requirement (i) therefore is fulfilled. In order 
to address items (ii) - (v) the collective behavior 
of an ensemble of ions has to be investigated. 

From part II) of the numerical simulations 
(5000 particles, short times) it is found that the 
population of ions in energy develops according to 
a maxwellian distribution function for a system of 
three degrees of freedom 



UE,t) 



exp 



2E^JHt) 



(19) 



where fa{E, t) has been normalized to n^, the den- 
sity of species a. 

By fitting this distribution model to the runs of 
part I) (500 particles, long times) the thermal en- 
ergy in time E*J^{t) = 1/2 mQ(w^''(i))^, and hence 
the heating rates, were extracted. The function 
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E^{t) for ^He and '*Hc is depicted in Figure 7. 
Clearly, '^He is accelerated faster than **He and 
this causes an enhancement at a given energy (say 
1 MeVamu"^). 
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Fig. 8. — Abundance ratio •^He/'*He of ions above 
1 MeV amu~^ vs. acceleration time. 

The ratio ^He/^He above 1 MeV amu~^ can 
now be calculated in dependence on acceleration 
time. We integrate 

oo 

nait)= J fa{E,t)dE , (20) 
1 MeV amu~^ 

where a = 3,4 corresponds to the according 
He isotope, and form the ratio (■^He/'*He)(t) = 
ns{t)/n4{t). It is found that the required abun- 
dance ratio of ^He/*He ~ 0.1 — 1 is reached 
in the time range of to ~ 0.63 — 0.89 s (see 
Fig. 8). Due to the higher acceleration rate of 
^He (about factor of 2-3 above ^He), the fast tail 
of the energy distribution function populates en- 
ergies above 1 MeV amu~^ faster than for ''He. 
This yields a ^He/^He » 1 for small times t <C toj 
^He/^Ho - 0.1 - 1 for to ~ 0.63 - 0.89 s and 
asymptotically approaches the coronal abundance 
ratio for times t ^ to- The derived time range 
of to corresponds well to the typical timescale of 
electron acceleration during the impulsive phases 
of flares of 0.1 — 1 s as manifest by the shortest 
hard X-ray peaks (Kiplinger et al. 1984). Since the 
EF wavefield is a direct consequence of the elec- 
tron acceleration process, the lifetime of the field 
is expected to be on the same order of magnitude. 
When electron acceleration ceases, the ^He/'^He 



ratio is frozen and no further ion acceleration oc- 
curs. 

The time integrated total number of ^He nu- 
clei being accelerated above 1.3 MeV amu"""^ in 
an impulsive flare is about 10"^^ particles (Reames 
et al. 1994). Assuming the flaring area to be of 
about 10^^ cm^ with a scale height of 10^ cm, 
the proton density to be = 5 x 10^° cm~^ 
and ^He/H = 5 x 10"^, there are 2.5 x 10^^ 
^He ions available in the flaring volume. This 
yields a percentage of about 4% of the pre-flaring 
^He population that has to be accelerated above 
1 MeV amu^^. Computing the value of the ■^He 
density n3(to) at time to and comparing it to the 
coronal abundance yields a percentage of about 

3.5 — 10.5%. Thus the number of accelerated 
ions at the required energies is sufficient to explain 
the observed particle fluxes. The model therefore 
also accounts for criterion (in). 
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Fig. 9. — Particle energy spectrum in MeV amu~ 
at time to = 0.76 s. ^He is harder than ^He , 
qualitatively reproducing the spectrum observed 
by Moebius et al. (1982). 

The ^He acceleration by EF waves also gener- 
ates ion energy spectra that qualitatively repro- 
duce the observed behavior. According to Moe- 
bius ct al. (1982) the observed spectrum of ^He 
is harder than the spectrum of "^He in the energy 
range of 0.4 — 4.0 MeV amu~^ which is a feature 
of all ^He rich periods in their study. The accord- 
ing particle densities in dependence on energy per 
nucleon resulting from the model presented herein 
are depicted in Figure 9, reproducing the observa- 
tions. 
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A turnover of the ^He spectrum at around ^ 
120 keV amu~^ energies can be seen in Figure 9. 
This reproduces very nicely the observations by 
Mason et al. (2000) who reported a turnover at 
around ~ 100 keV amu"^. It can be interpreted 
as the result of the whole ^He population being 
accelerated by EF waves. 

Therefore it has been established that the EFI 
is not only a viable acceleration mechanism for 
ions during solar flares, it also can account for 
the observed abundance enhancements of ^He over 
^He , accelerates enough ^He nuclei, and repro- 
duces the qualitative behavior of the particle en- 
ergy spectra. 

5. PROTONS AND HEAVIER IONS 

The energy density profile may suggest by the 
same arguments used to explain the enhancement 
of ^He above ^He, that H should be enhanced 
above ^He, which is not the case. This is not a real 
problem of the model because of the very differ- 
ent number densities of the according ion species. 
The energy in the wavefield available for proton 
acceleration is higher than the energy available for 
^He. Nevertheless, there are much more protons. 
Thus the available energy is distributed among a 
larger number of nuclei and the energy gained per 
nucleon is smaller for the protons. By the same 
argument, ^He should be further enhanced in com- 
parison to ''He. This effect was neglected herein 
since the scope of this paper is to show the ability 
of fractionation by the accelerator itself. 

When discussing the possible acceleration of the 
protons by EF waves, other effects have to be 
taken into account. Whereas, due to their low den- 
sities, ^He and "^He can be treated as test-particles, 
this view is not valid for protons and nonlinear ef- 
fects as the erosion of the spectral wave energy 
around the gyrofrequency of H have to be consid- 
ered. Such an analysis requires simulations that 
are beyond the scope of the presented work. 

"^He rich events usually are also characterized 
by enhanced Fe/0 ratios. Although the observed 
enhancement of Fe/0 is only about a factor of 
8 over the coronal ratio, it is still the strongest 
enrichment of the heavier ions. Other heavier ions 
as Ne, Mg, Si are enhanced by factors of 2-3 above 
O, whereas C, O, N are not enhanced. In this 
section the enrichment of Fe ions is addressed and 



briefly discussed in view of the model presented in 
this work. 

In a plasma at a temperature of lO'^ K, the most 
probable charge state of ^^Fe is 20-1- while "'^O is 
fully ionized. When resonantly cyclotron accel- 
erated, ^^O should exhibit the same acceleration 
rate as ^He while ^^Fe , following the argumenta- 
tion for enhanced ^He acceleration rate, is much 
slower heated. As can be seen from Figure 4c the 
energy available for ^^Fe is much smaller than for 
^^O . This behavior was confirmed by carrying out 
the same simulations for ^^Fe and ^^O as for ^He 
and ^He . Although the acceleration via EF waves 
by cyclotron resonance energizes ^^Fe ions up to 
energies of 1 MeV amu~^ on timescales of ~ 10 s, 
an enhancement of Fe/0 by the same mechanism 
as for ^He/^He presented above cannot be reached. 
It is in general impossible in the EFI model to ex- 
plain a simultaneous enhancement of ^He/^He and 
Fe/O when the spectral energy density exhibits a 
slope as depicted in Figure 4 c. The same energy is 
available for ^He and ^^O , but ^He has more and 
^^Fe has less energy available. The heating rate 
of ^^Fe will therefore always lie below the ^^O and 
*He rate, whereas ^He always lies above it in the 
present simple model. Other processes may enter 
the picture. Non-linear effects as the back reaction 
of the wavefield in response to energy loss to the 
ions can have a significant influence on the heat- 
ing rates. However, the analysis of such effects 
are beyond the scope of this paper and will not be 
addressed here. 

Similar problems arise for the enrichment of Ne, 
Mg, Si since their gyrofrequcncies also lie below 
the gyrofrequency of *He . TIk^ samc^ arguments 
as for ^^Fe therefore make an enhancement above 
^^O impossible in the present scenario. 

The observed heavy ion enrichment is indepen- 
dent of the degree of "^He enhancement, although 
heavier ion enrichments are generally associated 
with ^He (Mason et al. 2000). This is an indica- 
tion that it may not be the same mechanism that 
accounts for the acceleration of ^He and ^^Fe or 
other heavier ions. 

6. CONCLUSION 

A model for the enrichment of ^He during im- 
pulsive solar flares is presented. The accelera- 
tion of ^He and "^He can be understood as conse- 
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quence of Electron Firehose (EF) waves, resulting 
from an unstable anisotropic electron distribution 
with T^J > T^. The model does not need addi- 
tional sources of free energy or plasma properties 
than those resulting directly from the acceleration 
of the bulk of electrons. The essential result of 
the investigation is that EF waves, excited by an 
anisotropic electron distribution function, can ac- 
celerate ^He and *He via cyclotron resonance to 
MeV amu""'^ energies on timcscalcs of ^ 1 s. In 
support of this conclusion, the following specific 
results are found: 

1. The EF waves accelerate ^He ions via cyclo- 
tron resonance. The symmetry of EF waves in 
greatly enhances the efficiency of the acceleration 
process with respect to the standard cyclotron res- 
onant wave-particle interaction of a unidirectional 
propagating wavefield. The parallel force driving 
the ion out of resonance, and therefore limiting the 
acceleration time in the unidirectional case, can- 
cels out in the time average for the counterprop- 
agating wavefield. The total resonant interaction 
time therefore is strongly increased and the parti- 
cle can reach high energies. 

2. The linear growth of the EF waves self- 
consistently generates a spectral energy distribu- 
tion, causing enhanced acceleration of ^He with 
respect to ^He . It is found from test-particle sim- 
ulations that the heating rate of "^He exceeds the 
values for ''He by about a factor of 2-3. 

3. It is shown that this small enhancement in 
the heating rate of ^Hc above '^He already can ac- 
count for the observed enhancement in ■^He/''He 
from the coronal value of ~ 5 x 10~^ up to 0.1 — 1 
during impulsive solar flares. Due to the larger 
heating rates, the high energy tail of "^He populates 
energies above 1 MeV amu~^ faster. It reaches the 
required abundance ratio of '^Ile/^He for an accel- 
eration timescale of ~ 0.76 s. 

4. The fraction of ^He accelerated by the EF 
waves is large enough to account for the observed 
particle fluxes. Between 4 — 11% of the coronal 
^He population in the flare is accelerated above 
1 McV amu"^ . 

5. The acceleration model reproduces the qual- 
itative behavior of the '^Hc energy spectrum with 
respect to the ^Hc energy spectrum. As observed 
by Moebius et al. (1982) the ^He spectrum is gen- 
erally harder than the ^He spectrum. Moreover, 



a turnover in the ^He energy spectrum around 
~ 100 keV amu~^ is obtained, reproducing the ob- 
servations by Mason et al. (2000). 

EF waves are an inherent property of stochas- 
tic electron acceleration by transit-time damping 
of cascading fast-mode waves. Thus, their ability 
of accelerating ^He and ^He supports the scenario 
of stochastic acceleration of the bulk electrons in 
impulsive flares. Future work will have to extend 
these results to the oblique EF mode and investi- 
gate its role in the enhancement of heavy ions and 
in the flare acceleration process in general. 
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